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Abstract 

We study the stability of Bose-Einstein condensates with superfluid currents in a one-dimensional 
periodic potential. By using the Kronig-Penney model, the condensate and Bogoliubov bands are 
analytically calculated and the stability of condensates in a periodic potential is discussed. The 
Landau and dynamical instabilities occur in a Kronig-Penney potential when the quasimomentum 
of the condensate exceeds certain critical values as in a sinusoidal potential. It is found that 
the onsets of the Landau and dynamical instabilities coincide with the point where the perfect 
transmission of low energy excitations through each potential barrier is forbidden. The Landau 
instability is caused by the excitations with small q and the dynamical instability is caused by the 
excitations with q = n/a at their onsets, where q is the quasimomentum of excitation and a is the 
lattice constant. A swallow-tail energy loop appears at the edge of the first condensate band when 
the mean-field energy is sufficiently larger than the strength of the periodic potential. We find that 
the upper portion of the swallow-tail is always dynamically unstable, but the second Bogoliubov 
band has a phonon spectrum reflecting the positive effective mass. 
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I. INTRODUCTION 



Since its first achievement 1], Bose- Einstein condensates (BECs) in an optical lattice 
have been studied vigorously |2j. Various interesting phenomena have been observed in 



such systems, including the superfluid-Mott insulator transition 



flfl, 



gap solitons [5|, the Josephson effect 



formation o 

J 



brigh ; 
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Q], and breakdown of superfluidity 
121 ] . In particular, breakdown of superfluidity has attracted much attention. Recently, the 
stability of BECs in an optical lattice moving at constant velocity has been experimentally 
investigated and it was demonstrated that the Landau and dynamical instabilities occur 



when the velocity of the lattice potential exceeds certain critical values [12]. 



Swallow-tail energy loops h ave been t 
in a periodic potential 



14. 
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reoretically found in the band structure of BECs 



181 ] and attracted much interests. Mueller pointed 



out that a Bloch state in the upper portion of swallow-tails corresponds to a saddle point in 
the energy landscape and it is dynamically unstable In fact, Seaman et al. have shown 
the instability of a Bloch state in the upper portion of swallow-tails by numerically solving 
the time-dependent Gross-Pitaevskii (GP) equation [itJ. However, the stability of swallow- 
tails has not been fully investigated because a detail study of linear stability problem is still 
lacking. One needs to identify the unstable modes of excitations by solving the Bogoliubov 
equations. 

In this paper, in order to understand the superfluidity of BECs in an optical lattice and 
swallow-tail energy loops, we study BECs in a periodic potential by using the Kronig- Penney 
(KP) model. It is well-known in solid state physics that the KP model is useful for under- 

n 

standing the band structure of electrons [19|]. The KP model is also useful for the study of 
BECs in a periodic potential and has remarkable advantages. First, it allows one to cal- 
culate the whole band structure of condensate energy and excitation spectrum analytically. 
Another advantage is that the band structure of excitation spectrum can be calculated from 
the tunneling problem of Bogoliubov excitations through a single potential barrier. The 
tunneling problem of Bogoliubov excitations has been studied in several papers and it was 
found that low energy Bogoliubov excitations exhibit anomalous tunneling i.e. a potential 



barrier is transparent for them 20j, |2l|. In our previous work, the excitation spectrum of 



current-free condensates has been calculated using the KP model, and it has been shown 
that the —us tuuueuug is erueia, to the phoeon speet™ of exeitatious Q. 
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In the present paper, we extend the previous work to calculate the excitation spectrum 
of condensates with supercurrents. We solve the Bogoliubov equations using the solution of 
the tunneling problem of Bogoliubov excitations [23j, and discuss the stability of conden- 
sates with supercurrents. We will show that the Landau and dynamical instabilities occur 
when the quasimomentum of the condensate exceeds certain critical values. We identify the 
unstable modes causing the Landau and dynamical instabilities. We will also show that the 
Bloch state of the upper portion of swallow-tail is dynamically unstable, but the excitation 
spectrum has a phonon branch due to the positive effective mass. 

The outline of the paper is as follows. In Sec. II, we introduce our model and formulation 
based on the Gross-Pit aevskii mean-field theory. In Sec. Ill, we solve the time-independent 
GP equation and calculate the first condensate band, group velocity and effective mass. 
The condition for the presence of swallow-tails are discussed. In Sec. IV, we solve the 
Bogoliubov equations and calculate the band structure of the excitation spectrum. In Sec. 

V, we discuss the stability of BECs with superfluid currents in a periodic potential. In Sec. 

VI, we summarize our results. Some details of our calculation in Sec. Ill and IV are given 
in Appendix. 



II. MODEL AND FORMULATION 



We consider a BEC at the absolute zero of temperature confined in a combined potential of 
a box-shaped trap and a one-dimensional periodic potential. A box-shaped trap was realized 
in a recent experiment with a radial harmonic confinement and end caps The frequency 
of the radial harmonic potential is assumed to be much larger than the excitation energy 
in the axial direction so that the one-dimensional treatment of the problem is justified [22]. 
It is also assumed that the axial size of the system is so large that the effect of the edge of 
the system can be neglected. We consider condensates with supercurrents flowing through 
the periodic potential. This situation can be realized by an optical lattice consisting of 
two counter-propagating laser beams with a frequency difference and moving in a constant 
velocity 

HQ. 

The periodic potential is assumed to be an array of 5-function potential barriers, 

V(x) = V J2S(x-ja), (1) 
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FIG. 1: Schematic picture of a current-carrying condensate in a KP potential. The solid lines, gray 
bars and dashed arrow express the condensate density l^o^)! 2 ) potential barriers and supercurrent, 
respectively. 

where a is the lattice constant, and Vq is the potential strength. Note that Vq and a are 
independent parameters. Equation (0) is a special type of the KP potential. A schematic 
picture of the system is shown in Fig. ^ Although the periodic potential of optical lattices 
in experiments is sinusoidal, the KP model provides us with valuable insights. BECs in a 
KP potential have the same properties as those in a sinusoidal potential, such as the band 
structure Q and the excitation spectrum of a current-free condensate Q. As we 
discuss later, the Landau and dynamical instabilities of BECs occur in a KP potential as in 
a sinusoidal potential. 

Our formulation of the problem is based on the Gross-Pitaevskii mean-field theory j^ill^. 
The condensate wave function ^f(x, t) obeys the time- dependent GP equation, 



ih 



dt 



2m dx 2 



+ V(x) + g\V\ 



(2) 



where m is the mass of the atoms and q is the mean-field coupling constant. The coupling 



constant in ID is given by g = |2flj |. where ao is the s-wave scattering length and a± 
is the harmonic oscillator length of the radial confinement. The static solution ^q(x) of Eq. 
(0) satisfies the time-independent GP equation 

h 2 d 2 



- m + V(x) + g\* (x) 



* (x) = 0, 



2m dx 2 

where \x is the chemical potential. The normalization condition is given by 

>(j+l)a 



(3) 



W a (x)rdx = N, 



c- 



(4) 



.i<-i 



Nq is the number of condensate atoms per site. The energy of the condensate per site is 
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given by 



E 



0'+l)a 



dx% 



a 



2m dx 2 2 



9ur, ,2 



^0- 



(5) 



Considering small fluctuations from ^q(x), the condensate wave function can be written as 



e ft" 



^o(x) + u(x)e * — v (x)*e ft 



(6) 



Substituting Eq. (0) into Eq. (J2J) and keeping only linear terms of the fluctuation, (u(x), v (x)) 
obeys the Bogoliubov equations. 



*2 





h 2 d 2 



fi + V (x) + 2g\^ \ 



(7) 



(8) 



2m dx 2 

e/h and (u(x),v(x)) describe the frequency and amplitude of the collective mode of the 
condensate. 

Asserting the normalization condition 



/<MK*>l'-N.)l') 



1. 



(9) 



27, m 



one can quantize the collective modes and regard them as elementary excitations 
We hereafter call e and (u(x),v (x)) as the energy and the wave function of the elementary 
excitation. Strictly speaking, elementary excitations with complex energies do not exist, 
because the normalization condition Eq. is not satisfied if e is not real [3, 13, Q. 
However, we use this term by implicitly meaning collective modes with complex frequencies. 
The stability of BECs can be studied by calculating the excitation energy. The appearance 



of excitations with negative energies signals the Landau instability. It reveals that the 
solution \l/o of the time-i 
in the energy landscape 



solution \l/o of the time-indepenc ent GP equation does not correspond to a local minimum 



24L I25L |27j. This means that the process of spontaneous creation 
of excitations can take place and the system becomes unstable j^Jj. On the other hand, 
the appearance of excitations with complex energies signals the dynamical instability which 
means exponential growth of the fluctuation in time. The detail of the difference between 
the Landau and dynamical instabilities has been discussed in Refs. 
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III. CONDENSATE BAND AND SWALLOW-TAIL 



In this section, we shall obtain the Bloch state solution of Eq. (jSJ) and calculate the 
energy of the condensate E, group velocity v g and effective mass m* as functions of the 
condensate quasimomentum K. Assuming that the lattice constant is much larger than the 
healing length, analytic expressions of these quantities are obtained. Hereafter, we call the 
band structure of the energy of the condensate E(K) as the condensate band. 



A. Condensate wave function 

The condensate wave function can be written as ^o(%) — y/noA(x)e lS ^ x \ where A(x) and 
S(x) mean the amplitude and phase of the condensate. A(x) is normalized by the density 
at the center of each site n = \^o((j + 1/2) a) | 2 . Thus, Eq. (JHJ) is reduced to a set of 
equations as (we set h = 1 from now on) 

1 ''~ A + 9l A ~ 3 + V(x)A- fiA + gn A 3 = 0, (10) 



2m dx 2 2m 

Equation (fTT|) is the equation of continuity. Q describes the superfTuid momentum and the 
supercurrent is given by They are conserved for the static solution and independent 
of x. Note that the effect of the 5-function potential barriers appears only through the 
boundary condition at x = ja. 

Throughout the present paper, we assume that the condensate sits in the first condensate 
band. In this case, the condensate wave function has no node in a single well as shown 
in Fig. 2. 

Since the amplitude of the condensate wave function in the first band has a maximum at 
the center of each well, A(x) satisfies the boundary condition 



dA 

dx 



= 0. (12) 

Cj+1/2)o 



From the definition of the central density n , A(x) also satisfies 



A[[j + -)a)=l. (13) 
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FIG. 2: Amplitude A{x) of the condensate wave function when Vq = 5gno^o, a = 15^q, and 
Q = 0. 09£q 1 . The dashed and solid lines represent the exact solution Eq. (fi""")) and the approximate 
solution Eq. (|2(Jj) . respectively. 




FIG. 3: Phase S(x) of the condensate wave function when Vq = 5gno£o, a = 15£o> and Q = 
O.O^q 1 . The dashed line, thick solid line and thin solid line represent the exact solution Eq. (|18|) . 
approximate solution Eq. ((2*2*1) . and asymptotic solution Qx + sgn(x)<£>/2, respectively. 



The boundary condition at x = ja can be obtained by integrating Eq. ()10|) from ja — to 
ja + as 



dA 2 



dx 



dA 2 



ja+0 



dx 



AmV Q A 2 (ja). 



(14) 



ja-0 



Multiplying Eq. (jTUJ) by j£ and using Eqs. (fT*2*|) and ((TJIJ), Eq. (fTU|) can be integrated as 



2/i 



1)A 2 - A 1 



(15) 



Here, we introduced the healing length £ = { m 9 n o) Integrating Eq. (|TB*|) again .34], the 
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solution of Eq. (fTUj) in the region (j — |) a < x < (j + 1) a is given by 

y//3+ - P-(\x - ja\ + x ) 



A(x) 



1-/3 _W 



/3+-/3- 



where 



#n 




4(Q£o) s 



(16) 



(17) 



/i and x are determined as functions of V , a and Q by solving Eq. (fT3j) and Eq. (JT3J). Once 
the amplitude A(x) is obtained, the phase S(x) can be calculated as 



x Q 

dx—r^. 

J" 



'A 2 ' 



(18) 



S(x) - S(ja) 

The exact solution of Eq. (JTBJ) and lfTHj) has been obtained in Refs. 0,0]. The exact solution 
in a unit cell for (a, V ,Q) = (15 £ , 5 gn ^ , 0.09 ^q 1 ) is shown in Figs. El and El with dashed 
lines. From the exact solution, /i and xq are given by (/i, xq) = (1.0040520 grio, 0.1364985 £o)- 
For further analytic calculations, we assume that the lattice constant is much larger 
than the healing length. The condensate wave function far from each potential barrier 
asymptotically approaches the wave function of a uniform condensate. Accordingly, the 
chemical potential is given by 3 

H = gn + e Q} (19) 

where eg = Substituting Eq. ([TH]) into Eq. (JTBJ), an approximate solution of Eq. (fTUj) in 
the region (j — ~) a < x < (j + |) a can be obtained as 

A{xf = 7 (x) 2 + (Q£o) 2 , (20) 

where 

\A - (Q€o) 2 (\x - ja\ + x ) 



j(x) = Vl - (Q£o) 2 tanh 



(21) 



Substituting Eq. (j2*U|) into Eq. (fTgj) . the phase S(x) in the region (j — |) a < x < (j + |) a 
is given by 



S(x) - S(ja) 



ja 



dx% 
A 2 



Q(x — ja) + sgn(x 



tan 



i ( 



-tan 



(22) 
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From Eqs. (J2Uj) and (j22J), the condensate wave function is expressed as 

M%) = V^e^-^W^) ( 7 (x) - sgn(x) iQ^o) , (23) 

where 

>/7(j'a) 2 + (W 



The wave function Eq. (|23|) takes the similar form with that of a gray soliton [36|,|37|. The 
only difference between them is the constant xq in Eq. (|25|) which depends on the potential 
strength Vo. Substituting Eq. into Eq. (fTIj) . the condition for obtaining xq as a function 
of Q and Vo is 

7(ja) 3 + ^7(J«) 2 - (1 - (Q£o) 2 ) 7W + ^W^) 2 = °' ( 25 ) 

The approximate solution of the amplitude Eq. (J20|) and phase Eq. (|22|) in a unit 
cell are shown in Figs. El and El with solid lines. Vo and Q are chosen as (Vo, Q) = 
(5 gno^o, 0.09 ^o" 1 ). One sees that the exact and approximate solutions show a good agree- 
ment when a 3> £o- From the approximate solution, /1 and Xq can be obtained as 
(/i, x ) = (1.0040500 gn , 0.1364977 £ ) which are very close to the values calculated from 
the exact solultion. We will also apply this approximation to the calculation of the excita- 
tion spectrum in Sec. IV. Note that Eqs. (flfi*!) and pO]) coincide with Eqs. (11) and (15) in 
Ref. when Q = and x ^ 0. 



B. Condensate band, group velocity and effective mass 

According to Refs. the phase difference <p between the condensates in the neigh- 

boring sites is defined by 



<p = Q 



(i+l/2)a 



dx 



'(j-l/2)a 

Substituting Eq. into Eq. one obtains 



1 

I2 



(26) 



(p= 2 



tan 1 



1 - (QZ0) 2 



tan 



1 ( l{jo) 



(27) 



if can be interpreted as a phase jump across the potential barrier as shown in Fig. El ip and 
K are easily related by the Bloch theorem. 
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ip 

FIG. 4: Superfluid momentum Q as a function of the phase difference ip. The thick black line, the 
gray line and the thin black line represent the superfluid momenta with (a, Vq) = (15 £ av , 5 <?n av £ av ), 
) and (15£ av , 15#n av 6w), respectively. 

Substituting Eq. ()23j) into the normalization condition Eq. (J3J), one obtains 

1 - ^Vl-(^o) 2 + ^7(j'a) = — , (28) 
a a n 

where n av = is the average density. The healing length for the average density n av is 
denned by £ av = (mgn av ) _ 2. 

Equations (j2SJ), (j2*7j) and (|2~Hj) yield Q, -y(ja) and as functions of a, V and <p. If Eqs. 
(l23j) . (}2Tj) and (J2~%)) are expanded into the power series of gn av £ av /Vb and £ av /a assuming 
Vb » #n av £av and a > £ av , one obtains 

Q ~ ^sin^, (29) 
70'°) - ^ a ^ av (l + cos^), (30) 



#n av #ri av aVo 8V(f 



gn av a a 2 aV 

Substituting Eqs. (|2HJ), P^j) . (jHOj) and (pTTJ) into Eq. (jHJ), the energy of the condensate is 



A^av iV c ^ av 2aV 1 ^ 8V C 2 



+ ^TT^(1 - cosv.) + sinV, (33) 

o 



Ncgrisy 2 3a a 2 aVb 
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FIG. 5: Chemical potential /i as a function of the quasimomentum K. The thick black line, the 
gray line and the thin black line represent the chemical potentials with (a, Vq) = (15 £ av , 5 <7n av £ av ), 
(15£av,7.5 

S^avsav 

) and (15f av , 15#n av £ av ), respectively. 

0.60, 1 




FIG. 6: Energy E of the condensate per site as a function of the quasimomentum K, namely the 
first condensate band. The thick black line, the gray line and the thin black line represent the first 
condensate bands with (a,V ) = (15£ av , 5gra av £ av ), (15£ av , 7.5gre av £ av ) and (15£ av , 155n av £ av ), 
respectively. 

We have rescaled all the variables by the averaged density n av instead of the center density 
n , because n av (n ) is independent (dependent) of the condensate quasimomentum K. The 
higher order terms of gn av £ av /Vo and £ av /a are not shown just for avoiding the complication 
of the equations. The higher order terms of the expansion are shown in Appendix, because 
they are necessary for the calculations in the following sections. 

In Fig. |5 Q is shown as a function of (p. In Eq. (J29|) and Fig. HJ we can clearly see the 
well-known Josephson relation between the supercurrent and phase difference in the limit of 
Vb > #n av £ 

av I ■ 
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Imposing the Bloch theorem 

ty (x + a) = V Q (x)e iKa (35) 

on the condensate wave function Eq. (|23|). the relation between K and ip is derived as 

Ka = Qa + cp. (36) 

Substituting Eq. (jHHj) into Eqs. (J3*T|) and (jHB^) . /i(i^) and E(K) can be calculated as shown in 
Figs. ElandEl A loop structure of n(K) and a swallow-tail of E(K) are present at the edge 
of the first Brillouin zone for relatively small Vo (thick black lines). The loop structure and 
swallow-tail become smaller as Vo increases and they disappear for a certain critical value of 
Vq (gray lines). The condition for the presence of swallow-tail (loop structure) is obtained 
from Eq. (}36|) as follows. A swallow-tail is present if there are two values of ip which satisfy 
Eq. (|36|) with Ka = tt, while a swallow-tail is not present if there is only one value of tp 
which satisfies Eq. (I36j) with Ka = n. This condition yields the threshold 

Vo = (37) 

which agrees with the result obtained by the numerical calculations in Ref. 
The first derivative of the condensate band gives the group velocity, 

dE 



s dK 

Assuming Vo ^> gn av £ av and substituting Eq. (}33|) into Eq. (}38j) . v g is given by 



(38) 



n ~ 91/ Sill< ^' ^ 
C s ZVq 

where c s = a/^^P 1 is the sound velocity for a uniform condensate. Substituting Eq. (j3T?j) into 
Eq. (jnU), v g (K) is obtained as shown in Fig. [7| When Vo is small, v g (K) has a reentrant 
structure reflecting the presence of the swallow-tail (thick black line). When the potential 
strength is larger than the critical value of Eq. (|37|). the reentrant structure is absent. 

In the limit of Vq ^> gn^a, the first condensate band and group velocity take the form 
of the tight-binding approximation 391 l( 



E £7L=o , gn^t, 



2 



N c gn av N c gn av 2aV c 



c s 2Vb 

12 



9=0 +^^(l-cos^a), (40) 
2a Vq 

sin Ka. (41) 



5 av 







n 

Ka 



FIG. 7: Group velocity v g as a function of the quasimomentum K. The thick black line, the 
gray line and the thin black line represent the group velocities with (a, Vb) = (15£ av , 5 <?n av £ av ), 
(15£ av ,7.5#n av £ a v) and (15£ av , 15g-n av £ av ), respectively. 
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FIG. 8: Group velocity as a function of the quasimomentum K. The thick black line, the 
gray line and the thin black line represent the group velocities with (a, Vq) = (15£ av , 5<7n av £ av ), 
(15£ av ,7.5srn av £ av ) and (15£ av , 15#n av £ av ), respectively. 

One can define another type of group velocity by 



Assuming V 3> g^ av £ av and substituting Eq. (|3TJl into Eq. (J4l?j) . is given by 




-1.0, 



o 



Ka 



271 




(43) 




(44) 
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Ka 

FIG. 9: Effective mass m* . The thick black line, the gray line and the thin black line represent 
the effective masses with (a,V ) = (15f av , 5gn av £ av ), (15£ av , 7.5gre av £ av ) and (15£ av , 15#n av <!; av ), 
respectively. 

Substituting Eq. (pffij) into Eq. (}4*3|) . v^(K) is obtained as shown in Fig. |H1 diverges at 
the edge of the swallow-tail (thick black line). 

It is important to remark about the physical meaning of v g and v^. v g can be regarded 
as a superfluid velocity of a condensate flowing through a periodic potential, because n av v g 
is equal to the supercurrent On the other hand, the velocity of Bogoliubov phonons 

propagating in the opposite (same) direction to the supercurrent reduces (increases) b y v,, 
and that the Landau instability occurs when w M exceeds the Bogoliubov sound velocity [la, 
Q, 4llQ]- We will discuss the detail of this issue in Sec. V. 

The effective mass is defined by the inverse of the second derivative of the condensate 
band as 

-i / n \ -i 



m " = {mo) = [jg ■ (45) 



Assuming Vo gn av ^ av and substituting Eq. into Eq. (JEJ), m* is given by 

m* 2V OK 



(46) 



m gn av cosip dip 

Substituting Eq. (jHEJ) into Eq. pfi|) . m*(K) is obtained as shown in Fig. El The sign of the 
effective mass changes when the group velocity takes its extremum. m*(K) also changes its 
sign at the edge of the swallow-tail, because ^ diverges there. m*{K) is always positive in 
the upper portion of the swallow-tail and negative near the edge of the swallow-tail in the 
lower portion. This is consistent with the result of the exact solution Eqs. |T6|) and |T8|) . 
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FIG. 10: Condensate quasimomentum Kyi which gives the maximum group velocity 
47T/3 
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FIG. 11: Condensate quasimomentum which gives the edge of the swallow-tail. 
Since the sound velocity in a periodic potential is given by Q] 



Cb 



n R 



<9/i 



(47) 



m* dn av ' 

the negative effective mass leads to the dynamical instability due to long-wavelength 
phonons 



J41l, and the formation of bright gap solitons follows after the dy- 
namical instability jsl Lf| . Condensates are expected to be stable for the long- wavelength 
perturbation in the upper portion of a swallow-tail, because of the positive effective mass. 
However, the upper portion of a swallow-tail corresponds to energy saddle points in the two- 
state approximation and is predicted to be unstable 16]. The dynamical instability of the 
upper portion of a swallow-tail has been shown by numerically solving the time-dependent 
GP equation [l^j]. We will calculate the excitation spectrum in the following sections and 
show that the spectrum has a phonon branch but also has a different branch causing the 
dynamical instability. 

From Fig. 9, it can be seen that the effective mass changes its sign twice for positive 
K. Since the sign of the effective mass is closely related to the stability of condensates, it 
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is important to calculate the quasimomentum at which the effective mass changes its sign. 
Km is determined by the condition of maximum group velocity as 

dv„ 



dK 



0. (48) 



Assuming V 3> <7n av £ av , the phase difference y9 M when K = K M can be calculated from Eqs. 
(gEJ) and (IA12j) in Appendix as 

#n av £ av gn av ^ v (#n av £ av ) 2 
C ° SVM = — aV 2V*-- (49) 

Substituting Eq. (f4T)J) into Eq. (jSEJ), Km is obtained as shown in Fig. [THl Kmci asymptotically 
approaches | in the limit of Vo 3> gn av a, because the group velocity takes the form of Eq. 

SD. 

The effective mass changes its sign at the edge of the swallow-tail K = K^. K^ is 
determined by the condition, 

m* = 0. (50) 

Assuming V 3> gn av £ av , the phase difference (fs when K = K E is obtained by Eqs. (jBUjl and 
flA14|) in Appendix, 

cos^e = + —rj 2 ^T72 — + - T" 51 

gn w a V gn w a 2 2V Z a 2 

Substituting Eq. (J5T|) into Eq. (J3HJ), /Te is obtained as shown in Fig. [TTJ The size of the 
swallow-tail monotonically decreases until the potential strength reaches the critical value 

T/ gria, v a 

IV. BOGOLIUBOV EXCITATION SPECTRUM 

In this section, we shall solve the Bogoliubov equations and calculate the excitation 
spectrum of a current-carrying condensate in a KP potential. We hereafter call the band 
structure of excitation spectrum as the Bogoliubov band. As is well known, the band structure 



46]. 



of a single particle in a periodic potential can be calculated from a single barrier problem 
In our previous work, the Bogoliubov band of a current-free condensate in a KP potential 
has been calculated from a single-barrier problem for Bogoliubov excitations (3]. We shall 
extend this work to the case of current-carrying condensates. 

16 



A. Single-barrier problem 



In a recent work, the scattering problem of Bogoliubov excitations of a current-carrying 
condensate has been studied, and it has been found that the tunneling properties, such 
as transmission coefficient and phase shift, strongly depend on the presence of a supercur- 

n n 

rent [23j]. Here, we briefly review the calculation in Ref. |23j, since it is necessary for the 
calculation of the Bogoliubov band. 

We assume a single potential barrier V(x) = Vq5(x). The solution of the time-independent 
GP equation is given by Eq. (J23j) . The solution of the Bogoliubov equations with Eq. (j2Hj) 
can be calculated from the solution for a gray soliton obtained in Ref. j^. Substituting Eq. 
()23|) into Eq. (J7J), one can obtain the four particular solutions 



U n (x) = A ne i{(kn+Q)x+s S n(x)%] 




(k n £o) 2 gn 



2s 



2e 



i-(Q^r-i(xy + 



7(x) — i sgn(x) 

(k n £,o) 3 gn 



gno 



+ 



4s 



(52) 



v n (x) 



Kiogno 



2s 




(k n ^ ) 2 gn 



2s 



j(x) + i sgn(x) 

(knCofgno 



gn 



+ 



4s 



(53) 



for 1 < n < 4. k„ is a solution of 



Qk n 



m 



(54) 



k tp 

for a given s, where e^ n = j 2 -. Equation ()54|) is the Bogoliubov spectrum in a uniform 
system when the condensate flows with a velocity — j35j. (ui(x),Vi(x)) and (u2(x) , V2(x)) 
describe the scattering components. {uz{x < ti),v%{x < 0)) and (1*4(2; > 0),t>4(x > 0)) de- 
scribe the localized components around the potential barrier. {u±(x < 0),%(x < 0)) and 



(u 3 (x > 0),v 3 (x > 0)) describe the components which diverge for \x\ 
malization constant A n is given by 



oo 



The nor- 



A., 



e ia 2 



n 



n 



1,3,4 
= 2 



(55) 
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Supercurrent * 



* x 



Supercurrent ** 



FIG. 12: Schematic picture of the solutions of the single-barrier problem. Scattering processes 
in which (a) the excitation comes from the left-hand side, and (b) the excitation comes from the 
right-hand side. The solid curves, solid arrows, dashed curves, and dashed arrows describe the con- 
densate wave functions, scattering components, localized components and condensate supercurrent, 
respectively. 



where 



2gn e + + 2(gn - 2e Q )e kn + e k 



n rn 



+1 



2e^/gno (gn + e kn + e - ^) 



2eJgn (gn + e kn + e - ^) 



(56) 



The normalization constants of the scattering components are determined to satisfy 
(ui(x), vi(x)) — > (%,%) and (u2(x),V2(x)) — > (uk 2 ,Vk 2 ) for x — »■ oo, where the ampli- 
tudes Uk and Vk are the well-known Bogoliubov transformation coefficients given by 



Uk 



ign + t k + e 



2(e 



Qk^ 



Qk 

m c i(Qx+sgn(a:)^) 



'' //,n "' < /V "' : + m_ -i(QaH-sgn(x)f) 



2 e 



Qk 



(57) 
(58) 



Two independent eigenf unctions of Eq. (JJJ) corresponding to two types of scattering 
process are obtained by omitting divergent components. One is the process in which an 
excitation comes from the left-hand side (i/r(x)), and the other from the right-hand side 
(ip r (x)). Schematic pictures of the solutions of the single-barrier problem are shown in Fig. 
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0.1 0.2 



0.3 



FIG. 13: Transmission coefficient T 1 as a function of e and <p for Vo = 5grao£o- T l s with ip = 0, 
</? = 7r/3, y = o?a — 7r/2, 9? = 2-7r/3 and <p = tt are shown. This figure is almost the same as Fig. 



2(a) of Ref. 



23j where Vb = 10g?io£o- 



fT^l Here we consider the former written as 




+ a l 




x < 0, 



x > 0. 



(59) 



The coefficients r , a , t 1 , and 6 1 are the amplitudes of the reflected, left localized, transmitted, 
and right localized components, respectively. They are functions of the excitation energy e, 
the potential strength Vq and the phase difference of the condensate in the background tp. 
The boundary conditions at x = are given by 



^(+0) = ^(-0), 



dip 1 
dx 



+o 



dtp 1 
dx 



+ 2mV Q il> l (0). 



(60) 
(61) 



-o 



Equation (|6ip can be obtained by integrating Eq. (J7J) from —0 to +0. The coefficients r 1 , 



t 1 and b l are determined from Eqs. (|60|) and (|61|1. The analytical expressions of the 
coefficients were obtained in Ref. [23] in the case of e <C gn and Vq ^> gn ^ . Using the 
relations between the coefficients in ip l and i\f 



r\e, -<p) = r r (e, <p), a l (e, -(f) = a r (e, <p), 
t(e,-<p)=f(e,<p), b\e,-<p) = V{e,<p), 



(62) 



one can also calculate 
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0.1 0.2 

£ /(gfh) 



0.3 



FIG. 14: Phase shift 5 l of t 1 function of e and 99 for Vq = 5grao£o- 1° with ip = 0, = 7r/3, 
= ~ vr/2, 99 = 2-7r/3 and p = it are shown. 

In Figs. EUandEl the transmission coefficient T 1 = |t'| 2 and the phase shift 5 l = arg(t 1 ) 
are shown as functions of the excitation energy e and the phase difference ip. As shown in 
Ref. [23], the transmission coefficient strongly depends on the phase difference as follows. 
When ip = 0, T 1 approaches unity and 5 1 approaches zero as e approaches zero. This means 
that the potential barrier is transparent for the low-energy excitations. This behavior is 



called as anomalous tunneling 



spa: 

E 



2CJ, |2l(. As ip increases, the peak width of T decreases. The 



peak width becomes zero and the sign of 5 l changes when ip = <pa 



~ |. The anomalous 



tunneling is absent when (p = tp\. As ip increases further from the peak width increases 
and the anomalous tunneling appears again in the region ipA < P < ^ ■ 



The superfluid momentum Q takes its maximum value when p = ipA 
Pa is determined by the condition 



23j. Accordingly, 



dQ 

dcp 



0. 



(63) 



When Vq 3> gno^o, p>A can be calculated from Eqs. and (|A6|) (see appendix) as 



COS(^A 



2 

av 



(ffWav&A 

2V* 



(64) 



V aV 

Substituting Eq. (|54*j) into Eq. (jSHJ), Ka corresponding ip& can be calculated. Comparing 
Eqs. (|54*jl and (jUJ) , one can see that |.?GJ, at which the anomalous tunneling is absent, is 
always smaller than \Ku\, which gives the maximum supercurrent [4^. As we will discuss 
in Sec. V, K = Ka is the onset of the Landau instability in a KP potential. 
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B. Bogoliubov band 



Since the Bogoliubov equations are linear differential equations, a general solution of the 
equations can be described as a linear combination of independent solutions with the same 
energy. A general solution of the Bogoliubov equations with a single potential barrier can 
be written as a linear combination of ij) l {x) and ij) x {x) 



(65) 



where a and j3 are arbitrary constants. 

The solution of the Bogoliubov equations with the periodic potential Eq. can be 
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calculated by imposing the Bloch theorem to Eq. (J63j) as the current-free case in Ref. 
However, the non-diagonal element of Eq. (JJJ) does not have the periodicity of a when 
the condensate has a supercurrent, because ^q(x) satisfies the Bloch theorem Eq. (|35|). 
Introducing &o(x) and (ub(x),vb(x)) by 



\Kx 



V (x) = %{x)e 



(u(x),v(x)) = (u B (x)e iKx ,v B (x)e- [Kx ) 



(66) 
(67) 



the Bogoliubov equations reduce to 

H K -g$ {x) 

*2 



g$ {x 



-H 



K 




U B (X 



v B (x 



(68) 



H 



K 



1 ^ - K d \ T l \ n I* 12 

- i + e K - n + V (x) + 2g\%\ 2 



2m dx 2 m dx 
The non-diagonal element has the periodicity of a in Eq. 
be applied to iPb{x) = (ub(x), Vb(x)) as 



i/jb(x + a) = ipB(x)e 



iqa 



dipt 



dx 



dipi 



dx 



iqa 



(69) 

and the Bloch theorem can 

(70) 
(71) 



where q is the quasimomentum of the excitation. Solving Eqs. (fTUj) and (ITTj) at x 
one obtains 



2 ' 



exp 



. (—kx + k 2 )a 



+ (tf - rV r )exp 



. (&i - k 2 )a 



t exp 



-2q + ki + k 2 )a 



+ fexp 



. (2q - ki - k 2 )a 



(72) 
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The Bogoliubov band e(q) can be calculated by solving Eq. (J72J). Since we are assuming 
a £av, the localized components around the potential barriers, which decay exponentially 
and vanish at \x\ = |, do not appear explicitly in Eq. (J72l . 

If e is real, one can prove from the Bogoliubov equations that the Wronskian defined by 



Wt/SW) =u ] *^u i - u { ^-u j * +v j * - ^-v j * 
dx dx dx dx 



(73) 



is independent of x, where ipi and ip % are the solutions with the same energy. By evaluating 
W ip l ) and W(ip T *,ip T ), one obtains 



-h(\uk 2 


| 2 + 


\v k2 


\ 2 )-Q 


ki(\u kl \ 
h{\u kl \ 


2 +l 
2 + l 


v ki\ 


2 ) + Q 
2 ) + Q 


-h(\u k2 


l 2 + 


\ v k 2 


\ 2 )-Q 



J|2 



„r|2 



+ \t 



l|2 



.r|2 



(74) 
(75) 



which express the conservation law of energy flux 
obtains 



. By evaluating W{ip r *, ), one also 



{h (K| 2 + \v k2 ?) + Q) r l f* = (h (\u kl \ 2 + K| 2 ) + Q) r T *t l . 



(76) 



Substituting Eqs. (f7i|). (J75|) and (|75jl into Eq. (|72~j). the relation between e and g can be 
simplified as 



cos 


(ki~k 2 ) 
2 


1 2 









COS 



(2g -h- k 2 )a ^ -5 l + 5 V 



(77) 



where 



|t|e i5 , f 



\t e 



i5 r 



(78) 



It should be recalled that Eqs. (J72|l and (|77ji are valid only when a ^> £ av . Equation (|77|) is 
not valid when e is complex, because e is assumed to be real in the derivation. We have to 
solve Eq. (|72|) when complex exists. 

Solving Eq. (|77jl . the Bogoliubov band can be calculated. The Bogoliubov band 
for a = 15 £ av , Vo = 5 gn av £ av and if = is shown in Fig. lT5T a). Since is real and 
positive, the condensate with superfluid flow is stable when K = 4r. In contrast to the case 



of a current-free condensate in Ref. 



22|, e(q) is asymmetric with respect to q — 0. This is 



because the existence of the supercurrent breaks the left-right symmetry of the system. 
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FIG. 15: Bogoliubov band of a current-carrying condensate in a KP potential for a — 15<^ av , 
Vq = 55n av ^av arid K = The band structure is shown in (a) and the left-hand side of Eq. ((77|) 
as a function of energy is shown in (b). 



0.600 




Ka Km Kg. 



FIG. 16: Magnification of the swallow-tail energy loop in the first condensate band for (a, Vq) = 
(15£ av ,5#n av £ av ) shown in Fig. El 

In Fig. I15f b). the left-hand side of Eq. ()77|) as a function of e is shown. There exists no 
solution of Eq. (|77|) when the absolute value of the left-hand side exceeds unity, because the 
absolute value of the right-hand side is always less than unity. The energy regions for no 
solutions correspond to the band gaps. They are expressed as the shaded regions in Fig. [T^l 

V. STABILITY OF SUPERFLUID FLOW 

Since the phonons around q = are expected to be closely related to the stability, we 
shall calculate the excitation spectrum in the limit q —>■ 0. If one expands Eq. (|T2*|) assuming 
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q <C -, the linear dispersion of the excitation spectrum is obtained as 

(c b +<)|g|, q > 0, ^ 
(cb - Vn)\q\, g < 0, 

where Cb and are given by Eqs. (jUj) and (H2J), respectively. Equation (|75|) has been 
derived previously by using the hydrodynamic analysis the systematic g-expansion of 
the Bogoliubov equations |41j and the tight-binding model 
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42j. The change in the 



slope of the linear dispersion by is analogous to the change of the phonon dispersion in 
a uniform system due to the Galilean transformation. 

In Fig. E3 a magnification of a swallow-tail energy loop in the first condensate band is 
shown, and K&, Km and K^ are indicated by allows. We first consider a condensate with 
a positive group velocity v g in the lower portion of the swallow-tail (the region shown by 
the black solid line in Fig. [15)1 . The Landau (energetical) instability occurs if there exists 
any excitation with negative energy. One can clearly see from Eq. (f?H|) that the Landau 
instability occurs if > c s . We confirmed the fact that the onset of the Landau instability 
is due to excitations with small q by numerically solving Eq. (J72J1 . Therefore, the onset of 
the Landau instability is given by the condition 

v» = c s . (80) 

Assuming V 3> gn av ^ av and solving Eq. (jSDjl. one obtains 

gn av £ av . gn av £l v (gn av £ av ) 2 

coslfL - " + ~w 2^-' (81) 

where v^Lan is the phase difference corresponding to the critical quasimomentum for the 
Landau instability K Lan . From Eqs. (jo^j) and (J8~T|) . one sees that K Lan = namely the 
onset of the Landau instability coincides with the point at which the anomalous tunneling 
is absent. 

Dynamical instability occurs if there exists any excitation with complex energy. From Eqs. 
(I47|) and (|79J1 . one can see that excitations around q = cause the dynamical instability if the 
effective mass is negative, namely K > Km- On the other hand, it can be expected that the 
dynamical instability sets in at K < K M due to excitations around q = -, because crossing 
of phonon and anti-phonon branches, which is crucial to the appearance of excitations with 
complex energies, occurs initially at q = ^ j3, 13, Q . We calculate the critical value of K 
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FIG. 17: First Bogoliubov band with K < Ka for (a, Vq) = (15£ av , 5gn stv ( t£iv ). 

for the dynamical instability caused by excitations around q — -. The first Bogoliubov band 
is analytically expressed by assuming K ~ Ka- For avoiding the complication of equations, 
we introduce a variable b defined as 

b = — — — (cos<£>a — cosy?) , (82) 

av 

where —2 < b < 2. b is a function of i^, for example, 6(-^a) — and b(Kyi) = 2. Assuming 
Vo (^^av^av and substituting Eq. (|H2~J) into Eq. (|72^) . one obtains 

^^(-W sin2(9a) - W (?))' (83) 

Expanding Eq. (|H3^1 around g = — , it reduces to 



r2 



(tt - qa + v 73 ^) . (84) 



One can see from Eq. ()84)1 that the dynamical instability due to excitations with q — - 
starts at K = Ka- Therefore, both the Landau and dynamical instabilities start to occur 
at K = Ka- In the case of condensates in a sinusoidal potential, the onset of the Landau 
instability asymptotically approaches that of the dynamical instability in the limit of gn& v ^> 
-Er, where = is the recoil energy Since our assumption of sufficiently large 
lattice constant compared to the healing length corresponds to the condition gn^ 3> 
our result is consistent with that of condensates in a sinusoidal potential. 

In Figs.ElandEH the first Bogoliubov bands in the regions K < Ka and Ka < K < Kyi 
are shown, respectively. In these figures, one can clearly see that the Landau and dynamical 
instabilities start simultaneously at K = Ka- It is also clear that condensates in the region 
of Ka < K < K M are always unstable due to both Landau and dynamical instabilities. 
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FIG. 18: First Bogoliubov band with < K < Km for (a, Vq) = (15 £ av , 5gn sv £ iSCV ). (a) The real 
part of the first Bogoliubov band, (b) The imaginary part of the first Bogoliubov band. 

As K increases and approaches the slope of the linear dispersion reduces. The slope 
becomes zero at K = (Fig. IT7j) . When K exceeds excitations with negative energies 
appear (Fig. ITHT a)). At the same time, the imaginary part of the excitation energy grows 
around q — - (Fig. flST b)). For K > K M , all the excitations in the first Bogoliubov band 
have complex energies. 

We consider a condensate with a positive group velocity v g around the edge of the swallow- 
tail K = Kft, including the upper portion of it. In Figs. El and 1201 the first and second 
Bogoliubov bands around K = are shown, respectively. Since there always exist ex- 
citations with complex energies in the first Bogoliubov band (Fig. IT^|) . the upper portion 
of the swallow-tail is dynamically unstable. This resul t ag rees with that of the numerical 
simulation by the time-dependent GP equation in Ref. |l2(. On the other hand, we find in 
Fig. EDI that there also exists a phonon spectrum in the upper portion, reflecting the positive 
effective mass. As K in the lower portion of the swallow-tail approaches K^, the bottom 
of the second Bogoliubov band approaches the origin. The second Bogoliubov band turns 
into a gapless dispersion when K reaches K^, and it has a phonon spectrum in the upper 
portion. 
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FIG. 19: First Bogoliubov band around the edge of the swallow-tail for (a, Vq) = (15, 5). (a) The 
real part of the first Bogoliubov band, (b) The imaginary part of the first Bogoliubov band. 
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FIG. 20: Second Bogoliubov band around the edge of the swallow-tail for (a, Vq) = (15, 5). 



Thus, the stability of the condensate in a KP potential has been revealed in all the regions 
of the first condensate band. To conclude this section, we show the stability phase diagram 
in Fig. EU 



VI. CONCLUSION 



In summary, we have studied the stability and excitations of Bose-Einstein condensates 
with superfluid current in a Kronig-Penney potential based on the Gross-Pitaevskii mean- 

27 



K (quasi-momentum of condensate) 



FIG. 21: Stability phase diagram of the condensate in a KP potential when (a, Vb) = 
(15 £av, 5 <7n av c; av ). The light shaded area and the dark shaded area correspond to the regions 
of the Landau instability and the dynamical instability, respectively. This figure is focused to the 
swallow-tail energy loop. 

field theory. Solving the time-independent Gross-Pitaevskii equation, analytic form of the 
condensate wave function in the first condensate band has been obtained. Using this solution, 
the condensate band structure, group velocity and the effective mass have been calculated. It 
has been found that the first condensate band has a swallow-tail energy loop if 2Vn < gn w a. 
This condition agrees with the result obtained by numerical calculations in Ref. 17]. 

Imposing the Bloch theorem on the solution of the Bogoliubov equations in the single- 
barrier problem, we have calculated the Bogoliubov band and investigated the stability of the 
condensate. We have found that the onsets of the Landau instability caused by excitations 
around q = and the dynamical instability caused by excitations around q = ^ coincide, 
and anomalous tunneling of low-energy phonons for each potential barrier is absent at this 
point. When the effective mass is negative, condensates are dynamically unstable due to 
the excitations in the first Bogoliubov band. It has been found that the second Bogoliubov 
band has a phonon spectrum in the upper portion of a swallow-tail, while the upper portion 
is always dynamically unstable due to the excitations in the first Bogoliubov band. 
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APPENDIX: ANALYTIC EXPRESSIONS OF PHYSICAL QUANTITIES 

In this appendix, we show the higher order terms of the physical quantities expanded 
into the power series of gn av £ av /Vo and £ av /a- For simplicity, we introduce the following 
notation, 

c = cosy?, (Al) 



s = smip, (A2) 



a = A (A4) 

Sav 



K = KU. (A5) 

The superfluid momentum Q, j(ja), chemical potential /x, energy of the condensate E, group 
velocities v g and v^, and effective mass m* are expressed as 

* ^, s / 1 c 2 -l + 2c -2-c + 3c 2 \ 

Q = QU l + - + ^ + — + = + ^ , A6 

2V V a V a 2 aV 2V 2 J 

, He/ 1 -1-c 1 -5 + 3c -l-4c + 3c 2 \ 

7 ja ~ — ^- 1 + - + — = + — + ^- + =r- , A7 

2y V a Vb 2a 2 2aVb 2V 2 / 

w ~, ,„ v / 1 1 + c 3 5 + 5cc + c 3 



#n av \aV 8\/ 2 5 2 \/ 4fi\/ 2 4^ 3 

4 l + 3c -7 + 5c + 8c 2 -2 - 3c + 3c 2 + 4c 3 \ /A . 

a 3 V Q a 2 V Q 2 AaVJ 8V£ J 



2 2 2 1 6 8 4 2 . . . 

/i ~1 + -+ + _ + + A9 

a a 2 aV a 6 a 2 V a 3 V a 2 V<? aV£ 
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~ E ~, , . / 1 1 + c 2 1 + cc + c 3 

£ = ~ #L= + 1 - c — =- + - + — =- + — ^- + — 

A^c^av V2aV 8V 2 a 2 V W 2 4V n 3 



4 2 + 3c -5 + 7c + 10c 2 -2 - 3c + 3c 2 + 4c 3 \ /A _, 

H — H ~ 1 = 1 = , (A10) 

aW a?V 2 6W 3 8V* / 

^ =0 " 2 + 3^ + ^~^ + ^~^ + 3^~^ + a 2 ^ 2 + 35K?' ( U) 

v R dE s / 4 c 8 -2 + 3c -2 - c + 3c 2 . 
c s dK 2V V a K) a aV 2V 2 1 



v. da s (\ c 3 5c -1 + 2c - 3c 2 + 4c 3 4 

— = —hr ~ — ( 1 =- + — H — H = h — 

c s dK Vo \a W a 2 2W W 2 a 3 

-2 + 6c -6 - 3c + 12c 2 1- 12c-3c 2 + 16c 3 \ I dK 
+ ^V~ + 25^? + 8^ )/W ( } 



m 
m 



I I _c i_ _2c_ i -l+2c 2 , 4c , -3-2c+6c 2 , l-8c-2c 2 +9c 3 ' 

V / 2Vo SVo ^ V,? ^ 5?V ^ 2aV ~ l ~ 4V n 3 



where 



(A14) 



dK 1 c c -l + 2c 2 c -2-c + 4c 2 1 - 8c - 2c 2 + 9c 3 lk . 

— — ~ - + — — + — — H = 1 — H ~ 1 = . (A15) 

d<p a 2V W 2V 2 a 2 V 2W 2 W 3 

Equations ()A6jl . ()A12j) and (|A14J) are necessary for the calculation of Eqs. (JSIJ), (f4T)|) and 
(jHIJ, respectively. 
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